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Abstract 

The density fluctuation spectrum in the solar wind reveals a Kolmogorov-like scaling with a 
spectral slope of —5/3 in wavenumber space. The energy transfer process in the magnetized solar 
wind, characterized typically by MHD turbulence, over extended length-scales remains an unre- 
solved paradox of modern turbulence theories, raising the question of how a compressible mag- 
netofluid exhibits a turbulent spectrum that is characteristic of an incompressible hydrodynamic 
fluid. To address these questions, we have undertaken three-dimensional time dependent numeri- 
cal simulations of a compressible magnetohydrodynamic fluid describing super- Alfvenic, supersonic 
and strongly magnetized plasma fluid. It is shown that a Kolmogorov-like density spectrum can 
develop by plasma motions that are dominated by Alfvenic cascades whereas compressive modes 
are dissipated. 
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I. INTRODUCTION 



It is a curious observation that electron density fluctu ations in the int erstellar 
medium (ISM) exhibit an omnidirectional Kolmogorov-like Kolmogorov 19411 ] power 
spectru m k~ 5 ^ 3 (or -11/3 sp e ctral index i n three dimensions) o ver a 4 to 5 decade 



range [Armstrong et al 198ll . iHigdon 19841 . lArmstrong et al 19901 ] . The solar wind 
plasma also possesses density fluctuations that exhibit a Kolmogorov-like k~ 5 ^ 3 spec- 
trum (Goldstein et al 1995; Matthaeus & Brown 1988; Zank & Matthaeus 1990; 
Montgomery et al 1987; Lithwick & Goldreich 2001; Padoan & Nordlund 1999; 
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The Kolmogorov-like 5/3 spectrum is observed in many fluid, space and as- 
trophysical plasmas as well. For instance, turbulent spectra in the ISM and 
galaxies are foun d to exhibit a Kolmogorov-like sc aling in wa ven umber space 
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describing a Kolmogorov-like 5/3 spectrum in hydrodynamic fluid and magnetoplasma 
turbulence is however not po ss ible to cite here. Some go o d reviews such as 



Lesieur 1990 
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various possible scenarios of turbulence in general. While all of these works show the 
possiblity of a self-consistent energy exchange between widely disparate length-scales in 
the presence of waves, nonlinear structures, anisotropy, driving forces etc., the physical 
processes leading to a Kolmogorov-like turbulent density fluctuation spectrum is not yet 
fully understood. A great deal of attention has thus focused on understandi ng the evolu- 
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Higdon (1984, 1986) interpreted the observed density fluctuations Armstrong et al 1981 ] 
to be two-dimensional isobaric entropy variations in which temperature and density gradi- 
ents are directed oppositely and both are orthogonal t o the local app roximately uniform 
magnetic field. Based on a pseudosound approximation Lighthill 19521 ] . Montgomery et al 
(1987) related density fluctuations and incompressible magnetohydrodynamics (MHD) ve- 
locity and magnetic field fluctuations. This approach, called a pseudosound approximation, 
assumes that density fluctuations are proportional to the pressure fluctuations through the 
square of sound speed. The density perturbations in their model are therefore slaved to 
the incompressible magnetic field and the velocity fluctuations. This hypothesis was fur- 
ther contrasted by Bayly et al. (1992) on the basis of their 2D compressible hydrodynamic 
simulations by demonstrating that a spectrum for density fluctuations can arise purely as a 
result of abondoning a barotropic equation of state without even requiring a magnetic field. 
Bayly's (1992) work ignores magnetic field effects from the outset and hence does not ex- 
plain the influence and possible correlation of magnetic field and corresponding magnetized 
waves on the density fluctuation spectrum. The latter has been investigated in the slow solar 
wind plasma by Spangler & Spitler (2004) who suggest that there exists a strong correlation 
between the density and magnetic field fluctuations. The pseudosound fluid description of 
compressibility, justifying the Montgomery et al. (1987) approach to the density-pressure 
relationship, was further extended by Matthaeus and Brown (1988) in the context of a com- 
pressible magnetofluid (MHD) plasma with a polytropic equation of state in the limit of a 
low plasma acoustic Mach number (Matthaeus and Brown, 1988). The theory, originally 
describing the generation of acoustic density fluctuations by incompressible hydrodynamics 
(Lighthill, 1952), is based on a generalization of Klainerman and Majda's work (Klainerman 
and Majda, 1981, 1982; Majda, 1984) and accounts for fluctuations associated with a low 
turbulent Mach number fluid, unlike purely incompressible MHD. Such a nontrivial finite 
departure from the incompressibility state is termed a 'nearly incompressible' fluid descrip- 
tion and is put forward to provide a possible explanation of the turbulent density variations 
that are observed to exhibit a Kolmogorov-like power spectrum in the solar wind plasma 
(Montgomery et al. 1987; Matthaeus & Brown 1988, Zank & Matthaeus 1990, 1993, Shaikh 
& Zank 2004a,b). In the context of the ISM, a comparative study of two-dimensional turbu- 
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lence of self-gravitating supersonic MHD, hydro dynamic and Burgers turbulence by Scalo 
et al (1998) suggests a power law form of density fluctuations that is close to -1.7. Kim 
& Ryu (2005) report that the slope of the density power spectra in hydrodynamic turbu- 
lence becomes gradually shallower as the rms Mach number increases and it tends towards a 
Kolmogorov-like slope when the rms (or turbulent) Mach number is unity. The high resoul- 
tion simulations of Euler turbulence by Kritsuk et al (2007) suggest that the inertial range 
velocity scaling in the strongly compressible regime (with a spectral index close to -1.95) 
deviates substantially from the incompressible Kolmogorov 5/3-law. Kida & Orszag (1990) 
report that it is only the rotational component of the velocity field in driven hydrodynamic 
compressible fluid turbulence that exhibits spectra very close to that of the incompressible 
case even for a large Mach number (close to unity). Lithwick & Goldreich (2001) show that 
density fluctuations in the ISM plasma are generated by entropy modes while Alfven waves 
lead to a K olmoKorov-like spectrum. This, how ever, is not generic to all Alfvenic Mach num- 
bers (Ma) [Passot fc Vazquez-Semadeni 2003| . While both slow and fast modes introduce 
density fluctuations at large Ma, only t he slow mode dominates the densit y-magnetic field 
anticorrelation at relatively small Ma Passot fc Vazquez-Semadeni 20031 ] . However, the 
slow mode is known to be strongly Landau damped in a collisionless plasma. The presence 
of a mean magnetic field introduces additional complexities in the energy cascade processes 
in MHD turbulence. For instance, the assumption of isotropy breaks down in the presence 
the mean or large scale magnetic field. Along the direction of this large scale magnetic field 
Alfven waves suppress the parallel cascade (Iroshnikov 1963; Kraichnan 1965). The result- 
ing inertial range MHD spectrum is therefore thought to be flattened from fc -5 / 3 to k~ 3 ^ 2 
(Iroshnikov 1963; Kraichnan 1965). Within the paradigm of incompressible MHD turbu- 
lence, Goldreich & Sridhar (1995) proposed that parallel (k\\) and perpendicular (k±) modes 
are correlated by k\\ oc k 2 / 3 and energy in the perpendicular modes follows a Kolmogorov-like 
spectrum when linear (along B ) and nonlinear (across B ) frequencies balance. By contrast, 
Shaikh & Zank (2007) argue that such balance is not obeyed identically by the entire inertial 
range modes, but only by a few modes that critically balance linear and nonlinear Alfvenic 
frequency. Boldyrev (2005) proposes a scale dependent anisotropic power law that appears 
to differ from the in-situ solar wind observations (Podesta et al 2008). While in-situ space- 
craft observations of magnetic field fluctuations in solar wind plasma are shown to follow 
the Kolmogorov-like k~ 5 ^ 3 spectrum (Montgomery et al 1987; Zhou et al 1990; Goldstein et 
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al 1995), the velocity field tends to show a close consistency with a k 3 / 2 spectrum (Podesta 
et al 2006, 2007). The latter is contrasted by Roberts (2007) that the magnetic field and 
velocity in the solar wind do not evolve in the same way with helocentric distance. Based 
on Voyager observations, Roberts (2007a,b) argues that velocity spectrum relaxes towards 
a likely asymptotic state through spectral steepening and acquires a spectral index of -5/3, 
finally mathching the magnetic field spectrum. Roberts (2007b) further argues that -3/2 is 
accidental and transient, and that the -5/3 slope is the eventual state of all the fluctuations. 
In a comprehensive review, Bruno and Carbone (2005) and Veltri (1980) describe that low 
frequency solar wind velocity fluctuations closely follow a Kolmogorov-like spectrum, and 
the intermediate region (between high and low frequency) do not allow us to distinguish 
between a Kolmogorov spectrum (-5/3) and a Kraichnan spectrum (-3/2). Bavassano et 
al (2005) describe that large scale velocity field fluctuations in the polar solar wind closely 
follow a Kolmogorov-like 5/3 spectrum. 

The discrepancy in the magnetic and velocity field spectra continues to be an unresolved 
issue. Since our simulations tend to favor a Kolmogorov-like k~ 5 ^ 3 velocity spectrum in which 
density spectrum (A; -5 / 3 ) closely follows the velocity spectrum through a passive convection, 
we support a k~ 5 ^ 3 Kolmogorov-like spectrum for the velocity field fluctuations in solar wind. 
What is notable in our work is the dissipation of the high frequency component due to the 
damping (described in Section 2) of compressible plasma motion that suppresses the small 
scale and high frequency compressive turbulent modes. The MHD plasma therefore evolves 
towards a nearly incompressible state and the density is convected passively by the velocity 
field to yield a k' 5 ^ 3 spectrum. 

Perhaps, the most striking point about the Kolmogorov-like k~ 5 ^ 3 spectrum is its ubiq- 
uitous persistence in fluids and plasmas regardless of whether they are (in)compressible, 
(un)magnetized, (an)isotropic and (un)driven. The observed Kolmogorov-like density spec- 
trum yields two paradoxes; (1) why does a compressible magnetized fluid behaves as though 
it were incompressible and umagnetized, and (2) Why do the density fluctuations, an appar- 
ently quintessential compressive characteristic of magnetized turbulence, yield a Kolmogorov 
power law spectrum characteristic of incompressible hydrodynamic turbulence? These ques- 
tions have to be answered if we are to address the origin of the 5/3 spectrum in MHD 
turbulence in general and the ISM density power law spectrum, in particular. In this pa- 
per, we address these issues within the context of fully 3D simulations of compressible, 
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anisotropic, driven Magnetohydrodynamics (MHD) turbulence to understand how and why 
a supersonic, super Alfvenic, and low plasma (3 (f3 the ratio of plasma pressure and magnetic 
pressure) compressible MHD fluid should exhibit a Kolmogorov-like wavenumber spectrum 
in density. We find a strong correlation between the intrinsic MHD waves (i.e. Alfven, 
fast & slow magnetoacoustic modes) and nonlinear inertial range turbulent cascades that 
suggests that nonlinear mode coupling interactions in compressible MHD turbulence tend 
to dissipate high frequency magnetoacoustic modes. Consquently, the inertial range cascade 
is governed predominantly by Alfvenic interactions that passively convect the density field 
to yield a Kolmogorov-like k~ 5 ^ 3 spectrum. 

In section 2, we describe the governing equations, and critical assumptions of our 3D 
magnetohydrodynamic simulations. Section 3 describes results from nonlinear fluid simu- 
lations of compressible, driven, anisotropic, homogeneous, turbulent magnetohydrodynamic 
(MHD) plasma. Our simulations demonstrate that density, magnetic and velocity fields in 
MHD turbulence follow an omnidirectional Kolmogorov-like k~ 5 ^ 3 turbulent spectrum. The 
physical arguments explaining the evolution of a k~ 5 ^ 3 spectrum are described in section 4 
that deal primarily with turbulent damping of non-solenoidal velocity field fluctuations. In 
section 5, we outline our results for anisotropic cascades that result from the presence of a 
mean magnetic field in MHD turbulence. Finally, section 6 summarizes our major results. 



II. MHD MODEL 

Our underlying magnetohydrodynamic (MHD) model assumes that characteristic fluctua- 
tions in the magnetofluid plasma are initially isotropic, homogeneous, thermally equilibrated 
and turbulent. A large scale constant magnetic field is present and drives anisotropic turbu- 
lent cascades in an initially isotropic spectral distribution in compressible MHD turbulence. 
The characteristic turbulent fluctuations in the plasma are assumed in our model to be much 
bigger than shocks or discontinuities. In other words, we ignore the influence of shock on 
turbulent fl uctuations. Our work incorporati ng the effect of shocks on turbulent spectra is 



initiated in 



Zank et al 2007 



Zank et al 20 06] . The boundary conditions are periodic, hence 



mode coupling interactions in the local spatial region are considered. 

The fluid model describing nonlinear turbulent processes in the magnetofluid plasma, in 
the presence of a background magnetic field, can be cast into plasma density (p p ), velocity 
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(Up), magnetic field (B), pressure (P p ) components according to the conservative form 

dF r 



dt 



+ V • Qp = Q, 
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Equations (1) are normalized by typical length £ and time to = @o/Va scales in our 
simulations such that V = £ V,d/dt = t d/dt,V p = U P /V A ,B = B/Va(4vtp ) 1/2 , P = 
P/p V^,e p = e p /poV^,p = p/po- The bars are removed from the normalized equations (1). 
V A = J B /(47rp ) 1 / 2 is the Alfven speed 

The rhs in the momentum equation denotes a forcing functions (fjif(r, t)) that essentially 
influences the plasma momentum at the larger length scale in our simulation model. With 
the help of this function, we drive energy in the large scale eddies to sustain the magnetized 
turbulent interactions. In the absence of forcing, the turbulence continues to decay freely. 
While the driving term modifies the momentum of plasma, we conserve density (since we 
neglect photoionization and recombination). The large-scale random driving of turbulence 
can correspond to external forces or instabilities for example fast and slow streams, merged 
interaction region etc in the solar wind, supernova explosions, stellar winds in the ISM, 
etc. The magnetic field evolution is governed by the usual induction equation, i.e. Eq. 
(CE]), and obeys the frozen-in-field theorem unless dissipative mechanism introduces small- 
scale damping. Note carefully that MHD plasma momentrum equation contains dissipative 
terms on the rhs. It is the term in the momentum equation (i.e. d(p p U p )/dt- ■ •) that is 
proportional to /iV 2 Up+?7V(V-Up). The latter (along with the other terms in the equation) 
is divided by the density p p field during the evolution to calculate the velocity field. 



III. NONLINEAR SIMULATIONS 



We have developed a fully three dimensional compressible MHD code to study the nonlin- 
ear mode coupling inte raction in the context of compressible MHD t urbulence. Details of our 



code are described in 



Shaikh & Zank 2006 



Shaikh fc Zank 20071 ] . In the simulations, all 



the fluctuations are initialized isotropically with random phases and amplitudes in Fourier 
space and an initial spectral shape close to k~ 2 (k is the Fourier mode, which is normalized 
to the characteristic turbulent length-scale /o)- We have carried out the simulations for both 
decaying and driven-dissipative cases with and without external magnetic field Bq. Since we 
are interested in a local region of the solar wind magnetofluid plasma, the computational do- 
main employs a normalized three-dimensional periodic box of volume 7r 3 . Other parameters 
are 7 = 5/3, = 0.1 - 2.0, M A = 1.0 - 2.0, M s = 1.0 - 2.0, 77 = fj, = 10" 14 - 10" 15 , k f < 15.0, 
where 7, M4, M s , i], \x and kf are respectively the ratio of specific heats, Alfven Mach number, 
sonic Mach number, viscosity, magnetic diffusion and energy injection modes. Our MHD 
model does not include (photo) ionization and radiation terms. Hence source or sink terms 
corresponding to the density fluctuations are not included in our MHD model. Additionally, 
the localized dissipation is effective in our simulations at the small-scales where it damps the 
plasma motions. The small scale dissipation in the local inter stellar medium or solar wind 



may result from e.g. radiative cooling or ion-neutral collisions [Spangler 1991| . Accordingly, 
the small-scale dissipation in our model corresponds to collisional or viscous effects and is 
associated with the small-scale damping that is responsible only for the cascade of large- 
scale energy into the smaller scales thereby producing a well-defined inertial range turbulent 
spectrum. By contrast, the large-scales and the inertial range turbulent fluctuations remain 
unaffected by direct dissipation of the smaller scales. 

The initial kinetic and magnetic energies are equi-partitioned between the velocity and 
the magnetic fields. The latter helps treat the transverse or shear Alfven and fast/slow 
magnetosonic waves on an equal footing, at least during the early phase of the simula- 
tions. Magnetoplasma turbulence evolves under the action of nonlinear interactions in 



that larger eddies transfer th e ir energy to sma l 



cording to Kolmogorov 1941 



Iroshnikov 1963 



er ones through a forward cascade. Ac- 



Kraichnan 19651 ] . the cascade of spectral 



energy is mediated by a local interaction amongst the neighboring Fourier modes that 
continues until the energy in the smallest turbulent eddies is dissipated due to the finite 



S 



Reynolds number. This leads to the damping of small scale motions as well. This re- 
sults in a net decay of the turbulent sonic Mach number M s associated with the large 
scale fluctuations. If turbulence is not driven at large scales, the turbulent sonic Mach 
number continues to decay from a supersonic (M s > 1) to a subsonic (M s < 1) regime 



Shaikh fc Zank 20061 ] . This indicates that nonlinear cascades predominantly cause super- 



sonic MHD plasma fluctuations to become subsonic. In our decaying turbulence simulation, 
the large-scale energy simply migrates towards the smaller scales by virtue of nonlinear 
cascades in the inertial range and is dissipated at the smallest turbulent length-scales. On 
the other hand, spectral transfer in driven turbulence follows a similar cascade process 
as in the decaying turbulence case. However, the inertial range spectrum in the latter 
is maintained by a large scale forcing at k < 5. The spectral transfer of turbulent en- 
ergy in the neighboring Fourier modes in globally isotropic and homogene ous hydrodynamic 



and magnetohydrodynamic turbulence is the widely accepted paradig m 



that leads to Kolmogorov-like energy spectra, while Iroshnikov 1963 



Kolmogorov 1941 1 
Kraichnan 1965| de- 



scribe turbulent spectra in the presence of a mean or local Bq. The most striking ef- 
fect, however, to emerge from the decay of the turbulent sonic Mach number is that the 
density fluctuations begin to scale quadratically with the subsonic turbulent Mach num- 
ber as soon as the compressive plasma enters the subsonic regime, i.e. 5p ~ O(M^) 

when M, < 1. This was demonstrated in our 3D simulations of a compressible MHD 

f 1 

plasma [Shaikh &: Zank 20061 ] . It signifies an essentially weak co mpressibility in the magne; 



toplasma, and is consistent with a nearly incompressible state 



Zank fc Matthaeus 1990 



Zank fc Matthaeus 1993 



Matthaeus fc Brown 1988 . 



Shaikh fc Zank 2006 ] . 



In the context of the magnetoplasma being nearly incompressible, the density fluctations 
exhibit a weak compressibility in the gas and are convected predominantly passively in the 
background incompressible fluid flow field. This hypothesis can be verified straightforwardly 
by investigating the density spectrum which should be slaved to the incompressible velocity 
spectrum. This is shown in Fig. (TjQ) which illustrates that the density fluctuations follow 
the velocity fluctuations in the inertial regime over the long time (several Alfven transit 
time) evolution of MHD turbulence. The evolution of compressible magnetoplasma from 
a(n) (initial) supersonic to a subsonic or nearly incompressible regime is gradual and it 
results in the density field following the velocity fluctuations. In the subsonic regime, com- 
pressibility weakens substantially so that density fluctuations are advected only passively. A 
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FIG. 1: (Left) Velocity fluctuations are dominated by shear Alfvenic motion and thus exhibit a 
Kolmogorov-like /c -5 / 3 spectrum. The middle curve shows the magnetic field spectrum. Density 
fluctuations are passively convected by the nearly incompressible shear Alfvenic motion and follow 
a similar spectrum in the inertial range. The numerical resolution in 3D is 512 3 . (Right) The 
evolution of Alfvenic (&a) and fast /slow magnetosonic (&a/sO modes demonstrates that the spectral 
cascades are dominated by Alfvenic modes. 

passively convected fluid exhibits a similar inertial range spectra as that of its background 
flow field McComb 19901 ] . Likewise, subsonic density fluctuations in our simulations exhibit 
a Kolmogorov-like k~ 5 ^ 3 spectrum similar to the background velocity fluctuations in the 
inertial range. This, we believe, provides a plausible explanation for the Kolmogorov-like 
density spectrum observed in MHD turbulence i.e. they are convected passively in a field of 
nearly incompressible velocity fluctuations and acquire identical spectral features [as shown 
in Fig. (CQ)]. The passive scalar evolution of the density fluctuations is associated essen- 
tially with incompressiblity and can be understood directly from the continuity equation as 
follows. Expressing the fluid continuity equation as (d t + U • V) In p — — V ■ U, where the 
rhs represents compressiblity of the velocity fluctuations, shows that the density field is ad- 
vected passively when the velocity field of the fluid is nearly incompressible with V • U ~ 0. 
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The theoretical basis illustrating the nonlinear damping of the non-solenoidal component of 
velocity field, i.e. V ■ U — ► 0, is described quantitatively in the following section. 

IV. TURBULENT EVOLUTION OF NON-SOLENOIDAL VELOCITY FIELD 

Understanding, how an initially non-solenoidal velocity field evolves towards a solenoidal 
field is important as it explains the evolution of a compressible MHD magnetoplasma from 
a supersonic to a subsonic or nearly incompressible state that yields a passively advected 
Kolmogorov-like k~ 5 ^ 3 density spectrum. Any given velocity field can be decomposed into 
solenoidal and non-solenoidal components or, equivalently, longitudinal and transverse com- 
ponents, which we shall here refer to as fast/slow magnetosonic and Alfven components, 
respectively. The simulations show that the amplitudes of fast/slow magnetosonic compo- 
nents described by the quantity kus decay more rapidly than the amplitudes of the Alfvenic 
components. Consequently, they will be inefficient in cascading the corresponding inertial 
range spectral energy. Hence nonlinear interactions, in the inertial range, are governed 
predominantly by non-dissipative Alfvenic modes (kjC) that survive collisional damping in 
compressible MHD turbulence. This is quantitatively demonstrated in Fig. (TTJ) [right panel] . 

The damping of a non-solenoidal velocity component, in part, explains the observed 
Kolmogorov-like k~ 5 ^ 3 in our simulations. To this end, it is essential to distinguish the 
Alfvenic and non-Alfvenic, i.e. corresponding to the compressional or due to slow and fast 
magnetosonic modes, contributions to the turbulent velocity fluctuations. To identify the 
distinctive role of Alfvenic and fast/slow (or compressional) MHD modes, we introduce diag- 
nostics that distinguish the energies corresponding to Alfvenic and slow/fast magnetosonic 
modes. Since the Alfvenic fluctuations are transverse, the propagation wave vector is orthog- 
onal to the velocity field fluctuations i.e. k _L U, and the average spectral energy contained 
in these (shear Alfvenic modes) fluctuations can be computed as 



The above relationship leads to a finite spectral contribution from the |k x Uk| characteris- 
tic turbulent Alfvenic modes. On the other hand, fast/slow magnetosonic modes propagate 
longitudinally along the velocity field fluctuations, i.e. k || U and thus carry a finite com- 
ponent of energy corresponding only to the zk • Uk part of the velocity field, which can be 
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FIG. 2: Turbulent energy associated with the characteristic inertial range modes corresponding to 
|k x Uk| and k || U fluctuations is shown. Our simulations show that the velocity fluctuations are 
dominated by shear Alfvenic motion whose contribution, corresponding to the curve represented by 
|{7(/c)curiu| 2 ! is more than an order larger than that corresponding to the fast/slow magnetosonic 
modes and is shown by |£7(fc)Divu| 2 - This result is consistent with Fig 1 (right panel) that depicts 
the mode coupling evolution of the two MHD modes. 

determined from the following relationship 
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The expression of kMS essentially describes the modal energy contained in the non-solenoidal 
component of the MHD turbulent modes. 

The quantative evolution of the characteristic modes corresponding to the Alfvenic kA 
and slow/fast compressional magnetosonic &ms modes is depicted in Fig. 1 (right panel). 
Although the modal energies in kA and k^is modes are identical initially, a disparity in 
the cascade rate develops, and the energy in longitudinal (or compressional) fluctuations 
associated with the non-solenoidal velocity field decays far more rapidly than the energy in 
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the Alfvenic modes. The Alfvenic modes, after a modest initial decay, sustain the energy 
cascade processes by actively transferring spectral power amongst various Fourier modes. 
By contrast, the fast/slow magnetosonic modes (Icms) progressively weaken and suppress 
their corresponding spectral contribution in the turbulent energy cascades. The difference 
in the cascades corresponding to fc^ and kj^s modes persists even at long times. The 
/cms mode represents collectively a dynamical evolution of small-scale fast plus slow mag- 
netosonic cascades. The physical implication, however, that emerges from Fig. 1 is that 
the fast/slow magnetosonic modes (k MS ) do not contribute efficiently to the spectral trans- 
fer process, and that the cascades are governed predominantly by non-dissipative Alfvenic 
modes that survive the collisional damping in compressible MHD turbulence. Correspond- 
ingly, the turbulent energy associated with the Alfvenic modes makes a dominant contri- 
bution to the velocity fluctuation spectrum when compared to the magnetosonic modes. 
We have made measurements of the turbulent energy in the Alfvenic and fast / slow mag- 
netosonic modes to quantify their respective contributions to the velocity field fluctuation 
spectrum. This is shown in Fig. (2). In Fig (2), |t/(&)c U riu| 2 = Hfc |k x Uk] 2 / D/c ^ 2 and 
|t^(^)Divu| 2 — J2k |k ' Uk| 2 / Sfc k 2 . Clearly, the energy contribution by Alfvenic modes (par- 
allel to B and orthogonal to the velocity field) is more than 10 times that of the fast/slow 
magnetosonic modes (parallel to the velocity field). This clarifies that it is the predomi- 
nance of Alfvenic modes (Fig (1) & (2)) in inertial range cascades that primarily lead to 
a Kolmogorov- like spectrum. Damping of Alfven waves is possible by ion- neutral collisions 



as pointed out by [Kulsrud fc Pearce 1 969. 



Ba 



[for more references, see Shaikh fc Zank 20081 ]] . In our present simulations, we do not in 



sara 19961 ] in the context of molecular clouds 



elude ion-neutral damping as we focus mainly on the solar wind plasma. This nonetheless 
suggests that because of the decay of the fast /slow magnetosonic modes in compressible 
MHD plasmas, supersonic turbulent motions become dominated by subsonic motions and 
the nonlinear interactions are sustained primarily by Alfvenic modes thereafter; the latter 
being incompressible. One of the implications of the turbulent damping of non-solenoidal 
velocity field (ik • Uk) in an MHD fluid is that compressible modes make negligible or no 
contribution to the inertial range energy cascade. The cascade is thus determined primarily 
by the incompressible Alfven modes that passively convect the density fluctuations. This 
point is further consistent with the MHD fluid continuity equation which in k space reads 
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FIG. 3: (Left) Evolution of k\\/k± in anisotropic magnetofluid turbulence shows that k± becomes 
progressively dominant. (Right) Anisotropic magnetic field spectra along and across the Bq are 
consistent with the left panel. The spectra are computed in the steady state (close to 11 — 13 Iq/vq) 
and are from the same simulation. The forcing mode spans the wavenumber band 3 < kf < 5. 



as follows. 



(Jfi + i E S ( k + k ')U k ■ k'J In p(k, t) ~ -ik • U k . 



The nonlinear mode coupling interations associated with the Dirac delta function 5 are finite 
only for those interactions that obey the Fourier diad k + k' = in spectral space. It follows 
from our simulations that the turbulent damping of the non-solenoidal velocity field on 
the rhs of the continuity equation makes an insignificant contribution to the inertial range 
energy cascade. The nonlinear mode coupling interactions in MHD turbulence are therefore 
dominated by convective transport that leads to a passive convection of density fluctuations. 
The density fluctuations subsequently follow the inertial range spectrum and are identical 
to the backgrou nd nearly incomp r essible velocity field and thu s have a Kolmogorov-like 



k 5//3 spectrum Kolmogorov 194ll . iMcComb 1990l . iLesieur 19901 ] in MHD turbulence. Our 



simulation results described in Fig. (1) are fully consistent with this scenario. 
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FIG. 4: Spectral distribution of the ratio \B(k±)\ 2 /\B(k\\)f 



V. ANISOTROPIC TURBULENT CASCADES 



We find that the presence of a large scale magnetic field B (along the indirection) 
introduces anisotropy in the distribution of energy in wavevector space such that the rms 
wavenumbers along (k\\) and across (k±) the mean magnetic field Bq show a discrepancy i.e. 
h\ 7^ k± [see Fig. Q]. We employ the following diagnostics to monitor the evolution of the 
rms wavenumbers k± and k\\ in time. The rms k± mode is determined by averaging over the 
entire turbulent spectrum weighted by k±, thus 



k i 



(k ± (t)) 



E fc |fc±fl(M)l 2 
\ £ fc |£(M)l 2 ■ 
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Here (• • •) represents an average over the entire Fourier spectrum, k± = ^Jk% + ky. Similarly, 
the evolution of the k\\ = k z (along the B direction) mode is determined by the following 
relation, 

fcn = <*«(«)> 



Z k \k«B(k,t)\i 



\ E k \B(k,t)\i ■ 

The wavenumbers k x , k y and k z are respectively along the x, y and z directions. It is clear 
from these expressions that the k± and k\\ modes exhibit isotropy when k± ~ k\\. Any 
deviation from this equality corresponds to spectral anisotropy. We follow the evolution of 
k± and k\\ in our simulations. We find a disparity in the magnetic field fluctuation spectrum 
along 



and across 



TOi)I 2 = £|£(^ii)I 2 ^j 

k: 



\B(k ± )\ 2 ^J2\ B ( k ^ k \\)\ 2dk \\ 

fell 



the mean magnetic field (see Fig. 3, right panel). Note that the discrete summations J2k t 
and J2k ± are carried over k\\ and k± modes respectively. The presence of the mean magnetic 
field inhibits turbulent cascades in the parallel (to B ) direction and hence the characteristic 
modes (k\\) along the mean magnetic field are suppressed, while the modes in the orthogonal 
direction remain unaffected. The evolution in Fig. (3) therefore shows that the initial 
isotropic ratio k\\/k± ~ 1 progressively evolves towards anisotropy k\\/k± < 1. 

The suppression of the fen mode is caused by the excitation of Alfven waves, which 
act to weaken spectral transfer along the direction of propagation. This can be under- 
stood as follows; We assume that the spectral transfer, mediated by propagating Alfven 
waves, can be described by a three wave interaction mechanism, for which the frequency 



and wavenumber reso nance criteria are, respectively, expressed by Shebalin et al. 1983 . 
Matthaeus et al 1998 ] 

±Co> 3 —UJ\ — LU 2 , 

and 

k 3 = ki + k 2 . 

The frequency-wavenumber resonance conditions indicate that two Alfven waves (u>i, ki) and 
(u>2, k 2 ) mutually interact and give rise to the third wave k 3 ). Such conditions could, in 
principle, hold for a set of infinite waves as the indices '1' and '2' are merely dummy indices. 
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With the help of the Alfven wave dispersion relation {uo = k • Bo) and the component of 
wavenumber matching relation along Bq, i.e. k^ = ku +k2,,, we infer that either the fci,, = 
or fc 2 || = 0- K follows from the frequency-wavenumber resonance conditions that either the 
ki„ = or fc 2 ,, = 0. Owing to the absence of one of the parallel (to B ) components of the 
modes, nonlinear mode coupling interaction becomes inefficient in transferring the inertial 
range spectral energy in the parallel direction. Hence there is very little cascading along the 
magnetic field direction. Thus, the parallel wavenumbers (k\\) appear to be suppressed and 
the spectral cascade mainly occurs in the perpendicular wavenumbers {k±). Consequently, 
the magnetic field spectrum along the mean B is depleted. This, we suggest, explains the 
wavenumber disparity k\\ ^ k± observed in our simulations [see Fig 3]. While the turbulent 
cascades are effected locally by the presence of the mean magnetic field, the 3D volume 
averaged inertial range spectra in our simulations continue to exhibit a power law close to 
that of Fig. flQ. 

The question of anisotropic inertial range cascades has long been debat ed in the context 
of MHD turbulence. It dates back to the seminal work of Iroshnikov Iroshnikov 1963 ] 
and Kraichnan (1965), who first pointed out that the presence of a large-scale or self- 
consistently generated magnetic field influences the spec tral power cascade mechanism in 



a complicated manner 



Iroshnikov 1963 



Kraichnan 19651 ] . Kraichnan (1965) addressed the 



interaction of magnetized turbulent eddies with Alfven waves, excited as a result of a mean 
magnetic field. Owing to the presence of waves in a MHD fluid, turbulent correlations 
between velocity and magnetic field and the corresponding energy transfer time are deter- 
mined primarily by r ~ (bok)" 1 , where &o is a typical amplitude of a local magnetic field 
and is dimensionally identical to that of the velocity field by virtue of Elsasser's symme- 
try (Kraichnan 1965). According to Kraichnan (1965), it is this time scale that leads to 
the modification in the energy transfer (or cascade) because of the wave-turbulent eddy 
interaction process and without this modification the energy cascade rates would be deter- 
mined by typical hydrodynam ic eddy interaction time. The assumption of isotropy in the 



Kolmogorov energy spect rum 



large scale magnetic field 



Goldreich fc Sridhar 1995 



Kolmogorov 19411] was later modified to include the mean or 



Kraichnan 196 



5 



Shebalin et al. 1983 



Higdon 1984 



Higdon 1986 . 



Matthaeus et al 19981 ] . in part, because the presence of waves in 



MHD turbulence can significantly influence the energy cascade dyna mics in the wavenum - 
ber space. An attractive analysis was presented by Shebalin et al Shebalin et al. 19831 ] . 
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who demonstrated that the presence of a mean (or dc) magnetic field introduces a spectral 
anisotropy in the MHD turbulent spectrum. The observed anisotropic cascade in their work 
was understood to be due to the presence of Alfven waves, which are excited and propagate 
along the mean magnetic field. The non-dispersive propagating Alfven waves in the presence 
of a mean magnetic field give rise to distinct energy cascade rates along and across the large- 
scale magnetic field, thereby leading to an asymmetry in the spectral transfer rates. This can 
be understood as follows; The inertial range energy cascade rates depend on wavenumber (or 
modes). Since the presence of the background magnetic field depletes the parallel mode (i.e. 
fcii or k z ) and leaves the perpendicular mode k± unaffected, the spectral energy transfer cor- 
responding to the k« mode is suppressed in the direction of the mean or background magnetic 
field. By contrast, the perpendicular transfer of spectral energy remains unaffected. This is 
quantified by Fig (3), right panel in which the spectrum of \B{k\\) | 2 is suppressed as compared 
to that of \B(k±)\ 2 . It is because of this disparity, the energy cascades along and across the 
background magnetic field are distinct. Such asymmetric cascades, in agreement with ar- 
guments based on the frequency-wavenumber r esonance conditions d escribed as above, pro - 
duce spectral anisotropy in MHD turbulence Shebalin et al. 19831 . iMatthaeus et al 1998 ]. 
In the context of the solar wind plasma, fast (> 500km/s) an d slow (< 400km/ s) streams 
are dominated respectively by k\\ and k± anisotropic modes Dasso et al 20051 ] . A more 
quantitative treatment of anisotropic cascade was put forward by Goldriech and Sridhar 
Goldreich &: Sridhar 19951 ] for MHD turbulence, suggesting that the Kolmogorov energy 
spectrum corresponds to wavenumbers perpendicular to the local magnetic field so that 
E(k±) ~ k ± (where k± is wavenumber perpendicular to the local magnetic field), whereas 
the parallel wave num ber scales as fcj| oc k 2 / 3 . Bo ldyrev (2005) proposes that a Goldreich- 
Sridhar like spectrum Goldreich fc Sridhar 1995 1 results from a weak large-scale magnetic 
field, while the limit of strong aniso tropy, that is, strong large-scale m agnetic field, corre- 
sponds to the Iroshnikov-Kraichnan Iroshnikov 19631 . iKraichnan 19651 ] scaling of the spec- 
trum. This suggests that there exists an asymmetry in the energy cascade mechanism, 
which is believed to be due primarily to the presence of large-scale local magnetic fields in 
the turbulent MHD flow. Our simulation results describin g the anisotropic MH D turbulent 
cascades, depicted in Fig. (3), are further consistent with [Shebalin et al. 1983 ] . 
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VI. SUMMARY 



Our work proposes a self-consistent physical paradigm for the development of the den- 
sity fluctuations spectrum in solar wind plasma in the context of compressible, driven- 
dissipative, anisotropic 3D MHD turbulence. Understanding turbulent cascades in the 
MHD plasma is critical to many astrophysical phenomena. These range from understand- 
ing the role of waves and nonlinear cascades in the evolution of the solar wind, struc- 
ture formation at the largest scales, cosmic ray scattering and energization by solar wind 
turbulence at the smallest scales and the heating of the solar wind (Scalo & Elmegreen 
2004, Elmegreen & Scalo 2004, Goldreich & Sridhar 1995, Zank 1999, Goldstein et al 1995, 
Scalo et al 1 998, Goldreich 2001) to problem s such as energy transfer acros s many scales 
in the ISM Rov et al 20081. lHaverkorn 2008 



Elmegreen 1999 



Dickev et al 2001 



Rvu et al 20081 . iRickett 20071 . 



Willett 2005 . 



Minter & Spangler 1996 1. 



We find from our 3D (decaying turbulence) simulations that a fc -5 / 3 density fluctuation 
spectrum emerges in fully developed compressible MHD turbulence from nonlinear mode 
coupling interactions that lead to the migration of spectral energy in the transverse (i.e. 
k _L U) Alfvenic fluctuations, while the longitudinal "compressional modes" corresponding 
to k || U fluctuations make an insignificant contribution to the spectral transfer of inertial 
range turbulent energy. The explanation, in part, resides with the evolutionary characteris- 
tics of the MHD plasma that governs the evolution of the non-solenoidal velocity field in the 
momentum field. It is the non-solenoidal component of plasma motions that describes the 
high frequency contribution corresponding to the acoustic time-scales in the modified pseu- 
dosound relationship (Montgomery et al 1987; Matthaeus et al 1988; Zank Sz Matthaeus 
1990; Zank & Matthaeus 1993). What is notable in our present work is we find a self- 
consistent evolution of a Kolmogorov-like density fluctuation spectrum in MHD turbulence 
that results primarily from turbulent damping of non-solenoidal modes that constitute fast 
and slow propagating magnetoacoustic compressional perturbations. These are essentially 
a higher frequency (compared with the Alfvenic waves) component that evolve on acoustic 
time-scales and can lead to a "pseudosound relationship" as identified in the nearly incom- 
pressible theory (Matthaeus et al 1988; Zank & Matthaeus 1990; Zank & Matthaeus 1993; 
Bayly et al; Shaikh & Zank 2004a,b,c, 2006, 2007). The most significant point to emerge 
from our simulation is the diminishing of the high frequency component that is related to the 



19 



damping of compressible plasma motion. This further leads to the dissipation of the small 
scale and high frequency compressive turbulent modes. Consequently, the MHD plasma re- 
laxes towards a nearly incompressible state where the density is convected passively by the 
velocity field and eventually develops a k~ 5 ^ 3 spectrum. This physical picture suggests that 
a nearly incompressible state develops naturally from a compressive MHD magnetoplasma 
in the solar wind. 

The higher resolution in-situ solar wind velocity field observations at 1AU show a typical 
spectral index in the range from 1.4 to 1.5 (Podesta et al 2007; Tessein et al. 2009). The 
latter clearly differs from our simulations in which the spectrum of density field is close to 
1.67. To this end, we should point out that the solar wind observations at 1 AU contradict 
the fundamental theme of our model which states that the density fluctuations are convected 
passively by the velocity field and thereby acquire similar spectral power-law. Owing thus 
to the differences between the observations and our simulations, our model describing the 
passive convection of the density field through the nearly incompressible velocity field may 
not hold near 1 AU. Consequently, the solar wind density fluctuations (~ k~ 5 ^ 3 ) at 1 AU 
cannot be correlated with the velocity field (~ k~ 3 ^ 2 ) through the passive-advection phe- 
nomenology. We believe that the deviation of the density and velocity fields, leading thereby 
to the obvious differences near 1 AU, may arise from a number of physical processes, as noted 
below. Firstly, the physical processes driving the velocity spectrum near 1 AU are different 
from those beyond 1 AU. Therefore, the evolutionary characteristics associated with the 
velocity spectrum near 1 AU (describing a 3/2-like spectrum) can certainly not be represen- 
tive of the distant outer heliospheric turbulence spectrum. Secondly, close to 1 AU, stream 
interactions, shear instabilities, compressional modes, and many other physical processes 
can plausibly alter the velocity field spectrum. It is unclear from the work of Podesta et al 
(2007) and Tessien et al (2009) whether any of these processes have a significant influence on 
the observed spectral indices. To further clarify this point, we have carried out more simu- 
lations to distinguish the evolution of kinetic energy, and incompressible, compressible, and 
total velocity spectra (not shown in this paper). We find that the compressive velocity field 
exhibits a flatter spectrum, while the incompressible velocity field follows a 5/3 spectrum. 
This could very well mean that the velocity field, if is, driven or dominated by compressive 
modes, has a spectrum may be a flatter than 5/3. Thirdly, the Alfvenic interactions, often 
invoked to explain the descrepnacy between the Kolmogorov-like (5/3) and Kraichnan-like 
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(3/2) spectra, are ascribed to Alfven waves in MHD turbulence (Iroshnikov 1963; Kraichnan 
1965; Shebalin et al 1983; Biskamp 2003). The latter inhibits energy cascades along the di- 
rection of propagation. The spectral transfer of inertial range turbulent energy is therefore 
suppressed along the mean B , whereas the perpendicular cascade is governed predominantly 
by the hydrodynamic-like interactions. Consequently, the energy spectrum is dominated by 
hydrodynamic-like processes which lead to a Kolmogorov-like (5/3) inertial range turbulent 
spectrum. Hence a Kolmogorov-like (5/3) spectrum emerges in MHD turbulence when the 
nonlinear interactions are dominated by hydrodynamic-like eddies. By contrast, magnetic 
field eddies governing the Alfvenic interactions lead to a Kraichnan-like (3/2) spectrum in 
MHD turbulence. Thus, the disparate time scales associated with Alfvenic and compressive 
modes may flatten out the velocity field but not the density field. Our model, thus relating 
the density to the velocity field spectrum, may not hold near 1 AU in the circumstances 
where the velocity fluctuations are driven by the processes as described above. 

The Kolmogorov-like k~ 5 ^ 3 spectrum in density, velocity and magnetic fields resulting 
from our simulation (see Fig 1) are fully consistent with those of Tilley & Pudritz (2007), 
Mac Low et al (1998, 1989), Biskamp (2003), Padoan & Nordlund (1999), Stone et al (1998), 
Goldstein et al (1995), Goldreich & Sridhar (1995) and with others that are described 
elsewhere in our paper. 
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